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The excitation-chain theory of the glass transition, proposed in an earlier publication, predicts 
diverging, super-Arrhenius relaxation times and, via a similarly diverging length scale, suggests a 
way of understanding the relations between dynamic and thermodynamic properties of glass-forming 
liquids. I argue here that critically large excitation chains play a role roughly analogous to that 
played by critical clusters in the droplet model of vapor condensation. The chains necessarily induce 
spatial heterogeneities in the equilibrium states of glassy systems; and these heterogeneities may 
be related to stretched-exponential relaxation. Unlike a first-order condensation point in a vapor, 
the glass transition is not a conventional phase transformation, and may not be a thermodynamic 
transition at all. 



In (see also @), I proposed that both the dynamics 
and thermodynamics of glass-forming liquids might be 
determined by chainlike excitations of intrinsically dis- 
ordered molecular configurations. Because these excita- 
tions appeared to pertain only to irreversible behavior - 
e.g. to anomalously slow, super-Arrhenius, viscous relax- 
ation - it seemed remarkable that they could also govern 
equilibrium properties; thus I speculated that this situ- 
ation might imply partial violation of ergodicity at tem- 
peratures above the glass transition. After rethinking 
this situation, I now believe that such a drastic departure 
from conventional statistical mechanics is unnecessary. 
Rather, I think that the excitation chains allow us to 
understand the transition between a glass-forming liquid 
and a glass as an unconventional kind of phase transfor- 
mation or, quite possibly, not as a thermodynamic phase 
transformation at all. 

The most successful descriptions of the entropic col- 
lapse of glass forming liquids near the Kauzmann tem- 
perature have been spin-glass theories in which random, 
infinitely long-range interactions between Ising-like spins 
simulate glassy disorder. Such systems do not look much 
like liquids with short-range interactions between the 
molecules. Moreover, although such models produce 
qualitatively realistic thermodynamic behavior, it has 
been difficult to use them to compute super-Arrhenius 
dynamics. Bouchaud and Biroli 3], in reviewing this 
situation, have made an especially thoughtful effort to 
find a connection between thermodynamics and dynam- 
ics via the mean-field theories. In my opinion, their anal- 
ysis reveals deep problems associated with any attempt 
to describe localized phenomena such as droplets or do- 
main boundaries in a theory where the range of interac- 
tions is as long as the largest length scales in the system. 
Thus, I think it is useful to explore a different approach 
based only on short ranged interactions between ordinary 
molecules. 

Chainlike (or stringlike) excitations appear to be ubiq- 
uitous in the neighborhood of jamming transitions in 
amorphous materials. For example, see the work of 
Glotzer et al. A more recent example that seems rel- 



evant to the glassy state itself, just below the glass tran- 
sition, is seen in Fig. 4 of Silbert et al. || , which shows 
a low-frequency vibrational mode of a granular material 
near an un-jamming transition. The largest displacement 
amplitudes lie predominantly along stringlike paths. As 
in [1| , I visualize an excitation chain in a glass as a string 
of molecular displacements producing the glassy equiva- 
lent of a "vacancy" at one end and an "interstitial" at 
the other. That is, an excitation chain is a momentary 
density fluctuation that most probably lasts for only a 
few molecular vibration periods, but may under some 
circumstances induce a long lasting change in the sys- 
tem, moving it from one inherent state to another. It 
is the latter possibility, where the vacancy and intersti- 
tial in effect dissociate from each other, that is relevant 
to alpha relaxation processes above the glass transition. 
These chainlike excitations also must be present in the 
glassy state below the transition, where they cannot en- 
able the system to explore different inherent states but 
must be present and contribute - perhaps only negligible 
amounts - to the thermodynamics because they are ther- 
mal fluctuations just like any other thermal excitation of 
the system. 

There is an interestingly imperfect analogy between 
the excitation-chain model of a glass and the droplet 
model of a vapor in the neighborhood of its condensation 
point. In the latter, one approximates the vapor as a non- 
interacting gas of liquid droplets. The partial pressure as- 
sociated with droplets consisting of n molecules, when n 
is large, is assumed to be po(n) oc exp(— fn/ksT) where 
/„ is the free energy of a size-n cluster: /„ « h n + cr n 2 / 3 . 
Here h is the chemical potential measured from the two- 
phase coexistence point and a is proportional to the sur- 
face energy. In the stable phase, i.e. for h > 0, pu(n) is 
an exponentially decreasing function of n for arbitrarily 
large n. Any cluster of molecules that reaches a sub- 
stantial size n is more likely to evaporate than to grow. 
On the other hand, in the metastable phase with h < 0, 
there is a critically large n beyond which po{n) begins to 
increase, meaning that droplets larger than this critical 
size are more likely to grow than to evaporate. This crit- 
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ical droplet is the transition state for nucleation; thus, 
historically, the droplet model has been used primarily 
for estimating condensation rates. Equilibrium thermo- 
dynamic properties (away from critical points) are deter- 
mined primarily by small clusters, i.e. low-order terms in 
the virial expansion, for which the droplet approximation 
makes little sense. There also are intrinsic uncertainties 
in computing /„ for realistic situations even for very large 
n; in particular, it is not clear how to decide which clus- 
ters of molecules are close enough together to be counted 
as droplets. Such uncertainties have largely been elim- 
inated in field-theoretic calculations by looking only in 
the immediate neighborhood of the transition state, i.e. 
at the "instanton" [fjl]]], and computing the flux through 
this point in function space to find the condensation rate. 
The instanton theory also predicts a branch cut in the an- 
alytic continuation of the free energy from the stable to 
the metastable phase. The discontinuity across the cut 
is proportional to the condensation rate; thus there is a 
deep connection between equilibrium and noncquilibrium 
behaviors in these systems. 

The excitation-chain model of a glass has much the 
same structure as the droplet model of a vapor. In 
I argued that the probability of thermal activation of 
an excitation chain consisting of N molecular steps and 
extending a distance R (in units of the average molec- 
ular spacing) is proportional to exp(— AG(N, R)/k B T), 
where AG(N, R) is a free energy in the sense that it in- 
cludes a sum over all the allowed configurations of the 
chain. AG(N, R) consists of several parts: 



AG(N, R) = AG C 
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The first term, AGoo, is the bare activation energy, i.e. 
the energy required to form the vacancy and the inter- 
stitial regardless of whether or not these defects are well 
separated from each other. The remaining terms on the 
right-hand side of Eq. JJ) describe the excess free energy 
of the chain. ATeo, is the bare activation energy of the N 
links of the chain. The average energy per link, eo, is the 
energy required to move a pair of molecules far enough 
away from each other to allow a third molecule to pass 
between them. Ei nt is an energy that makes it unfavor- 
able for the links of the chain to lie near each other. Here 
I have used Flory's approximation for the self-exclusion 
energy of a polymer chain: 



E mt (N,R) a k B T mt 
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The last term on the right-hand side of Eq.GJ, 
W(N,R), is a sum over chain configurations. In |l| I 
wrote W in the form: 



In W(N,R) a uN-^-^-R, 



(3) 



7(T) = 7o(T /T) 2 is the mean-square strength of the 
fluctuations in e /k B T. Here, T = e a /vk B is the char- 
acteristic temperature determined by the energy eo; and 
7o is an inverse localization length associated with the 
fact that the excitation chains exist in the highly dis- 
ordered environment characteristic of a molecular glass. 
Because lengths are measured in intermolecular spac- 
ings, and glassy disorder is on molecular scales, 70 must 
be of order unity near the glass temperature. (I omit 
the diffusion term, R 2 /2N, from the usual random-walk 
analysis because it is negligible for present purposes.) 
Note that both the Flory approximation in Eq.© and 
the configuration-counting approximation in Eq.© make 
sense only in the limit of large N and R. 

The next step in this analysis is to evaluate the for- 
mation probability for all chains of length N by find- 
ing the value of R = R* at which A G(N, R) is a 
minimum, and then approximating this probability by 
pxc(N) oc exp(-AG*(iV)/fc B T), where 

AG*(N) = AG(N,R*) = AGoo - Nk B v{T -T Q ) 
+4 (V6) 3/4 (k B Tint) 1 / 4 (l(T) k B T) 3 / 4 N 1 ' 2 . (4) 

Note the similarity between pxc(N) and the formula for 
Pr>(n) in the droplet model. For T < T , pxc(N) de- 
creases monotonically for all N. However, for T > T , 
AG*(iV) has a maximum, and pxc(N) a minimum, at a 
value of N, say N*, which I identified in Q as the length 
of the critically large chain that nucleates a long-lasting 
density fluctuation. If a chain fluctuates to a size larger 
than N*, it is highly likely to continue growing with- 
out bound, and thus to dissociate the vacancy-interstitial 
pair. The activation energy for this process is 



AG* (T) AG{N*,R* 



^ + ^(T)R*(n (5) 



where exp (v) is the number of choices that the suc- 
cessive links in the chain can make at each step, and 



k B T k B T 
where the length R*(T) is 

Equations (JSJ and © recover the Vogel-Fulcher law for 
the logarithm of the a relaxation time r Q , which diverges 
linearly at T . Thus the excitation-chain mechanism pro- 
vides a simple explanation for super-Arrhenius behav- 
ior in glasses, and also predicts a diverging length scale 
R*(T). 

My thesis is that the excitation chains play a role near 
the glass transition analogous to that played by large 
droplets near a condensation point. The equilibrium par- 
tition sum in the glass phase, below To, includes a dilute 
population of long excitation chains that, like the large 
droplets, contribute negligibly to the thermodynamics 
but are present nonetheless. As in the droplet model, 
the overwhelmingly most important contributions to the 
thermodynamics come from small, localized, thermal dis- 
placements of the molecules, which may sometimes re- 
semble short chains in a glass. We certainly do not want 
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to use the long-chain formulas to estimate the contribu- 
tions of the short chains. Just as in the droplet model 
where we cannot decide which molecular clusters are to 
be counted as droplets, there is no unambiguous defini- 
tion of a thermally excited chain; we do not know how 
large the molecular displacements along the chain must 
be in order for a sequence of such displacements to be 
counted as a "chain." The best we can do is guess that 
these displacements are of order the molecular spacing 
for critical chains of size N* , and then suppose that this 
guess will become more precise if we can construct the 
analog of an instanton theory for this system. 

The crucial point is that large, rare, excitation chains 
are intrinsically part of the equilibrium state of a glassy 
system, both below and above To, and that they play an 
especially important role above To where they nucleate 
transitions between inherent states. With this under- 
standing, we may deduce that the thermal equilibrium 
state of a glass-forming liquid is necessarily heteroge- 
neous on a length scale of order R* (T) . A homogeneous 
glassy region much larger than R*(T), according to the 
formula for pxc(N), would contain a substantial popu- 
lation of chains longer than N* and thus would be un- 
stable against changes in size and shape. Conversely, an 
isolated domain smaller than R* could not support exci- 
tation chains and therefore would be immobile. Presum- 
ably, the molecules in the interiors of these domains are 
frozen into energetically favorable glassy configurations 
determined by their frustrated, short-range interactions 
- just as in the glassy state below T . The boundaries 
between these domains must then consist of more disor- 
dered material. In the language I have used here, they 
contain a high density of vacancies and interstitials. Pic- 
tures of such a domain structure have appeared in nu- 
merical simulations. For example, see Figs. 6b and 7 in 
I 

Note that the excitation chains cannot cross the dis- 
ordered boundaries, where the vacancies and interstitials 
at the ends of the chains can rccombine with preexisting 
interstitials and vacancies. These recombinations cause 
the domain boundaries to wander slowly on the time scale 
r a . As a result of this wandering, there is no breaking of 
ergodicity above T . The system samples all of its state 
space, but does so increasingly slowly as T decreases to- 
ward To. The glass- forming liquid freezes continuously 
as T decreases; the frozen domains become larger and 
the unfrozen molecules in the boundaries between the 
domains become a vanishingly small fraction of the sys- 
tem. Over long times at fixed T, any individual molecule 
is sometimes in a frozen state (inside a domain, locked 
into a cage by its neighbors, and contributing to the en- 
tropy as if it were in the glassy state below To ) , and is 
sometimes unfrozen (in a boundary region, contributing 
to the entropy as if it were more nearly in a liquid state). 
At any given time, the unfrozen fraction of the molecules 
is proportional to the surface-to-volume ratio of the do- 



mains, 3/R*(T). Thus the configurational entropy per 
molecule, i. e. the excess over the residual entropy in the 
frozen glass, in units of fcs, is 



Sc(T) 



3 sp 
R*(T) 



6\ 1/2 s ^(T-T ) 



it) (7oT T irlt ) 1/2 ' 
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where so is the configurational entropy per unfrozen 
molecule, multiplied by the molecular thickness of the 
inter-domain boundary, which I guess is independent of 
the size of the domains. 

Equation (JJJ implies that the excess configurational 
entropy vanishes linearly in T at the Kauzmann temper- 
ature Tk, and that Tk = Tq. The combination of Eqs.JSJ) 
and yields 



AG* (T) 7T7QS0 



k B T 



■XT) ' 



(8) 



which is the Adam-Gibbs formula. 9] 

Definitions of the glass temperature T g generally have 
the form AG*{T g )/kBT g = X g , where X g is a number 
of order 30, chosen roughly to represent the observable 
limits of long relaxation times or high viscosities. Then, 
because AG*(T g ) diverges near T , T g = T , the fragility 
to is LLCl 



_ d ( AG*(T) 

m = - T df 



T=T„ 
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Returning to the thermodynamic formula, Eq.J7Jl, we find 
that the jump in the specific heat at T g = Tq is 



T 



dsc 
dT 



T=T 



6 s 7q m 



(10) 



Here we recover the conjectured proportionality between 
Ac p and to, but with two material-specific parameters 
that might account for the observed scatter in the exper- 
imental data, so should scale with the "bead" number 
of the molecules, which ordinarily is factored out in ob- 
taining the linear relation in Ea. IjlOfl . [lll | Because it is 
basically a geometrical quantity, 70 may be roughly a 
constant, of order unity. This analysis also implies that 



R*(T g 



7o ' 



(11) 



Thus the critical length scale R* at the glass tempera- 
ture is predicted to be independent of the fragility, in 
qualitative agreement with results shown by Berthier et 
al. |l2j : but those authors report a substantially smaller 
length scale. 

The droplet-model analogy suggests that some kind 
of transformation is occurring at To, but it is not clear 
whether this is a true phase transition or, if so, what 
kind of transition it might be. Dissociation of a vacancy- 
interstitial pair does not nucleate a qualitatively new 
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state of the whole system, as does the appearance of 
a critically large droplet in a supersaturated vapor; it 
is simply a mechanism by which the glass takes steps 
in exploring its configuration space and thus starts to 
behave like a liquid. As discussed in there is a 
super- Arrhenius range of temperatures, Tq < T < Ta, 
within which the excitation-chain mechanism is opera- 
tive and the a relaxation time has approximately the 
Vogel-Fulcher temperature dependence. Ta is the tem- 
perature at which N* and R* become small, i.e. where 
excitation chains are no longer needed to stabilize den- 
sity fluctuations. Above Ta, relaxation is Arrhenius - the 
activation energy is temperature independent - and the 
system behaves like an ordinary, unjammed liquid, albeit 
like one where molecular rearrangements still must over- 
come activation barriers. The transition from the glass to 
the liquid occurs continuously across the super- Arrhenius 
region, with no latent heat, nor any change in symmetry. 

The excitation-chain picture does not require that To 
or Ta be sharply defined temperatures. The divergence 
of AG*(T) at T = T in Eq. © may be an artifact 
of the approximations made in deriving Eq.(|3J). (See Q 
for a discussion of the uncertainties at low T and large 
N, where the discreteness of the glassy molecular struc- 
ture must become increasingly relevant.) The actual re- 
laxation rates may remain finite but unobservably small 
down to T = 0, in which case there would be no thermo- 
dynamic phase transition at all. In principle, ergodicity 
would be restored at all temperatures, but thermody- 
namic equilibrium well below T g would be experimentally 
inaccessible. Such a possibility would be consistent with 
the recent results of Donev et al. 0], who find that a 
binary mixture of hard disks cannot exhibit a transition 
to an ideal glassy state of zero entropy, although that 
system does behave in a glasslike manner at sufficiently 
high densities. 

Similarly, we know that the large-A^ approximations 
used in Eqs. (J2J and are invalid for small N and 
cannot produce the smooth transition between super- 
Arrhcnius and Arrhenius behaviors that is seen exper- 
imentally at Ta- In all likelihood, Ta has no special 
thermodynamic significance, although it does mark the 
approximate upper bound of the temperature range in 
which heterogeneities on the scale R*(T) are present. 

The spatial and temporal heterogeneities implied by 
this theoretical picture make it likely that a closer look 
at relaxation processes will reveal stretched exponential 
decays of correlations. A detailed discussion of this topic 
is beyond the scope of the present report, but one illus- 
tration may be useful. (The following analysis has some 
features in common with the trapping models that have 
been discussed, for example, in references [lil fisl Hfij|V 
For T just above To, any individual molecule spends most 
of its time immobilized within a frozen domain. Occa- 
sionally, a boundary region wanders through its position 



and allows it to diffuse a substantial distance before be- 
ing frozen again. If the molecule is frozen in a domain 
of size R, the probability of its remaining in place for 
a time t (in units of r a ) is proportional to exp(— t/R 2 ). 
Then, if the distribution of domain sizes is Gaussian with 
a width R*, ~ exp[— (R/R*) 2 ), the most likely value of 
R scales like i 1 / 4 , and the relaxation function becomes 
cxp(— 2i 1 / 2 /i?*). If this crude estimate is valid, the lim- 
iting value of the stretching exponent would be 1/2. As 
T increases toward Ta, the domain size decreases and the 
exponent must return to unity in a smooth way that - like 
other features of the transition between super- Arrhenius 
and Arrhenius behavior - remains to be determined in 
this theory. 

In summary, I believe that the excitation-chain hy- 
pothesis may provide a framework - but not yet more 
than a framework - for constructing a theory of the glass 
transition based on realistic molecular models. Out- 
standing challenges include improved estimates of acti- 
vation rates for both long and short chains, a systematic 
statistical mechanical derivation of the entropy formula 
in Eq.lJZJ, and an analysis of how and where stretched- 
exponential relaxation might emerge from this model. 
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